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1.  INTRODUCTION 


The  external  propulsion  accelerator  was  proposed  by  Rom  (1990)  and  certain  operating  characteristics 
were  presented  by  Rom  and  Kivity  (1988)  and  Rom  and  Avital  (1992).  The  external  propulsion 
accelerator  was  proposed  for  develqxnent  based  on  the  successful  development  of  the  ram  accelerator  by 
Hertzberg  and  colleagues  at  the  University  of  Wa^ngton  (Hertzberg,  Bruckner,  and  Bogdanoff  1988). 
Several  issues  related  to  the  theory  and  operation  of  the  external  [Hopulsion  accelerator  are:  1)  the 
establishment  of  combustion  by  a  small  step  on  the  projectile  and  the  possibiliQr  of  flame  quenching  on 
the  projectile  afteibody,  2)  the  generation  of  raough  thrust  to  counter  the  drag  produced  by  the  step, 
3)  centering  of  the  projectile  within  the  accelerator  tube  at  launch  and  the  aeroballistic  stability  during 
flight  through  the  tube,  and  4)  potential  advantages  of  the  external  propulsion  accelerator  as  compared  to 
the  ram  accelerator.  This  report  addresses  the  first  two  issues  and  leaves  the  remaining  issue  for  further 
study. 

In  the  ram  accelerator,  a  projectile  and  obturator  are  injected  at  supersonic  velocity  into  a  stationary 
tube  filled  with  a  pressurized  mixture  of  hydrocarbon,  oxidizer,  and  inert  gases.  Flow  stagnation  on  the 
obturator  initiates  combustion  of  the  mixture.  A  system  of  shock  waves  generated  on  the  projectile  and 
reflected  from  the  accelerator  tube  wall,  in  conjunction  with  viscous  heating,  sustains  combustion.  The 
resulting  energy  release,  which  travels  with  the  projectile,  also  generates  high  pressures  that  impart  thrust 
to  the  projectile.  The  U.S.  Army  Research  Laboratory  (ARL)  ram  accelerator  system  uses  a  120-mm  (bore 
diameter)  tube  that  is  modeled  after  the  38-mm  system  at  the  University  of  Washington  (Hertzberg, 
Brackner,  and  Bogdanoff  1988).  Numerical  simulation  conducted  at  ARL  using  finite-rate  chemical 
kinetics  has  been  used  to  investigate  this  flow  field.  Numerical  results  are  used  to  visualize  the  flow  field, 
predict  effects  of  variation  in  system  parameters,  and  predict  projectile  in-bore  velocity  (Nusca  1994). 
Ram  acceleration  (30-mm  and  90-mm  bore)  has  also  been  demonstrated  at  the  Institute  of  St  Louis  (ISL) 
in  France  (Giraud,  Legendre,  and  Simon  1991;  Smeets  et  al.  1994). 

Unlike  the  ram  accelerator,  thrust  for  the  external  propulsion  accelerator  is  generated  by  combustion 
that  is  initiated  by  a  normal  shock  wave  generated  on  the  projectile  body  using  a  ramp  or  fcHward-facing 
step.  Since  the  projectile  is  subcaliber,  the  process  is  independent  of  the  tube  wall  (except  that  the  tube 
contains  die  pressurized  mixture).  There  are  potentially  difficult  design  issues  such  as  the  ability  to  initiate 
and  stabilize  the  combustion  front  on  the  projectile  in  such  a  way  that  positive  thrust  can  be  generated. 
The  step  that  initiates  combustion  also  produces  a  drag  component.  Previous  studies  (Rom  and  Kivity 
1988;  Rom  and  Avital  1992)  indicated  that  if  a  relatively  small  forward-facing  step  on  the  projectile 
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shoulder  (i.e.,  the  foiebody-afterbody  junction)  is  able  to  induce  a  detonation  wave  and  a  sustained 
combustion  front,  then  reasonable  thrust  levels  may  be  achieved. 

The  present  investigation  examines  the  ability  of  a  small  forward-facing  step  to  ignite  the  premixed 
fiiel/oxidizer  mixture  flowing  at  hypersonic  speeds  over  the  projectile  (due  to  projectile  movement).  The 
fiowfield  and  combustion  front  generated  by  the  interaction  of  the  bow  shock  from  the  projectile  nose  and 
the  detonation  wave  generated  by  the  step  are  investigated.  This  interaction  results  in  the  confinement 
of  a  high  {Hcssure  and  temperature  region  below  a  contact  surface  (see  Figure  1).  Ihis  region,  located 
between  the  projectile  body  and  the  contact  surface,  acts  like  the  combustion  chamber  of  an  engine. 
However,  in  this  case,  the  contact  surface  acts  as  the  "engine  cowling"  without  a  structural  boundary.  The 
high  temperature  and  pressure  combustion  products  are  then  expanded  over  the  afteibody  of  the  projectile 
to  produce  thmst.  The  flowfield,  with  chemical  reactions,  is  studied  in  detail  using  the  numerical  solution 
of  the  Reynolds  Averaged  Navier-Stokes  (RANS)  equations,  including  the  equations  for  chemical  species 
and  finite-rate  kinetics. 


Figure  1.  Schematic  of  projectile  and  reacting  flowfield. 

The  present  investigation  uses  the  solution  of  the  RANS  equations  for  flows  with  chemical  processes 
ad£q>ted  for  the  ram  accelerator  configuration  at  ARL  by  Nusca  (1991, 1993, 1994)  and  Kraczynski  and 
Nusca  (1992)  and  is  siunmarized  in  this  repoit.  In  these  calculations,  the  same  mixture  of  CH4,  O2.  and 
N2  selected  at  ARL  for  the  ram  accelerator  firings  is  used.  The  chemical  reaction  mechanism  for  this 
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mixture,  used  in  the  calculations  for  the  ARL  ram  accelerator,  has  yielded  good  results  when  compared 
with  experimental  data  (Nusca  1994).  Although  the  combustion  conditions  in  the  case  of  the  external 
propulsion  accelerator  are  different  from  those  in  the  ram  accelerator  (i.e.,  scale  of  the  projectile  and 
method  of  ignition),  it  has  been  assumed  in  this  report  that  by  using  the  same  mixture  with  the  same 
kinetics  model,  plausible  results  can  be  achieved  for  the  external  propulsion  accelerator  as  welL  This 
report  presents  calculations  for  the  effects  of  step  height  and  projectile  Mach  number,  using  several 
projectile  geometries,  on  the  combustion  process  and  resulting  thrust  These  calculations  are  also  used 
to  study  the  flow  ahead  of  the  forward-facing  step,  including  the  formadtm  of  a  combustion/detonatitxi 
front 

2.  REACTING  FLOW  MODEL 

Computational  fluid  dynamics  (CFD)  flow  simulations  for  the  ram  accelerator  projectile  were 
performed  at  ARL  using  the  Rockwell  Science  Center  USA-RG  (Unified  Solution  Algorithm  Real  Gas) 
code  (Chakravarthy  et  al.  1985;  Palaniswamy  and  Chakravarthy  1989;  Palaniswamy,  Ota,  and 
Chakravarthy  1991).  Nusca  (1991,  1993,  1994)  and  Kruczynski  and  Nusca  (1992)  used  this  code 
successfully  at  ARL  for  simulation  of  a  full-bore  12()-mm  ram  accelerator  projectile.  This  CFD  code 
solves  die  full,  3-D,  unsteady  RANS  equations,  including  equations  for  chemical  kinetics  (finite-rate  and 
equilibrium).  These  partial  differential  equations  are  cast  in  conservation  fonn  and  converted  to  algebraic 
equations  using  an  upwind  finite-volume  formulatioa  Solution  takes  place  on  a  mesh  of  nodes  distributed 
in  a  zonal  fashion  around  the  projectile  and  throughout  the  flowfield  such  that  sharp  geometric  comers 
and  other  details  are  accurately  represented.  Ibe  conservation  law  fonn  of  the  equations  assures  that  the 
end  states  of  regions  of  discontinuity  (shocks,  detonations,  deflagrations)  are  physically  correct  even  when 
smeared  over  a  few  computational  cells.  The  total  variation  diminishing  (TVD)  technique  is  employed 
to  discretize  inertia  terms  of  the  conservation  equations,  while  the  viscous  terms  are  evaluated  using  an 
unbiased  stencil.  Flux  computations  across  cell  boundaries  are  based  on  Roe’s  scheme  for  hyperbolic 
equations  (Roe  1981).  Third-order  spatial  accuracy  can  be  maintained  in  regions  of  the  flow  field  with 
continuous  variation  while  slope  limiting,  used  near  large  flow  gradients,  reduces  the  accuracy  locally  to 
avoid  spurious  oscillations. 

The  RANS  equations  for  2-D/axisymmetric  reacting  flow  (N  species  mixture)  are  written  in  the 
following  conservation  form  (Palaniswamy  and  Chakravarthy  1989)  for  dependent  variables  of  energy  (e), 
density  (p),  axial  and  radial  momentum  flux  (pu,  pv),  and  species  flux  (pa). 
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In  diese  equations,  Cj  and  0|  arc  the  mass  fraction  and  chemical  {n-oduction  term  for  the  i-th  species, 
respectively.  The  species  viscosity  (p)  and  themal  conductivity  (k)  are  referenced  to  ,  and 

using  Sutherland’s  law. 
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where  Tj, ,  Tq  ,  S„,  and  vary  with  species  (see  Stull  and  Prophet  1971).  The  mixture  viscosity  and 

P  K 

thermal  conductivity  are  determined  using  Wilke’s  law  (Wilke  1950)  denoting  f  as  p  or  k, 
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where  and  Mj  are  the  mole  fraction  and  molecular  weight  of  the  i-th  species,  respectively.  Pick’s  law 
is  used  to  relate  the  mixture  diffusivity  to  the  mixture  viscosity  through  the  Schmidt  number  Sc  = 
Pn/(pD).  The  specific  heat,  enthalpy,  and  Gibbs  free  energy  of  each  species  (per  mass)  are  givai  by  the 
following  fourth-order  temperature  (T)  polynomial  curve  fits  (Drummond,  Rogers,  and  Hussaini  1987) 
(AHg  is  the  heat  of  formation  and  Rj  is  the  specific  gas  constant  for  the  i-th  species). 

.5^  =  Aj  +  BjT  +  CjT^  +  DjT^  +  EiT'*  (9) 
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The  mixture  enthalpy,  total  energy  per  unit  volume,  and  ratio  of  specific  heats  are  given  by 
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Cp.  =  lEf-l  /%idT 


(15) 


A  covolume  equation  of  state  is  used  to  relate  pressure  (p)  to  temperature  (T). 


p(V-b)  =  nR„T,  b  =  ]C?.inibi, 


(16) 


where  V  is  the  voliune,  b  is  the  covolume  (Hirshfelder,  Curtiss,  and  Bird  1954),  n  is  the  number  of  moles 
of  gas,  and  R„  is  the  universal  gas  constant. 
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The  system  of  Navier-Stokes  equations  in  equation  1  is  valid  for  the  laminar  flow  of  a  viscous 
Newtonian  fluid.  In  reality,  the  flow  will  remain  laminar  up  to  a  certain  critical  value  of  the  Reynolds 
number,  pUL/^,  where  U  and  L  are  representative  values  of  the  velocity  and  length.  Above  this  critical 
value  the  flow  becomes  tuibuloit  and  is  characterized  by  the  aiq)earance  of  fluctuations  in  all  variables 
(U,  p,  p,  T,  etc.)  around  mean  values.  These  fluctuations  are  statistical  in  nature  and  cannot  be  described 
in  a  deterministic  way.  In  this  investigatitm,  the  Reynolds-averaging  approximation  is  used  along  widi 
algebraic  tuibulence  models  that  relate  fluid  viscosity  (p)  to  other  variables.  See  Nusca  (1993)  for  further 
details. 

Hydrocarbon  mixtures  such  as  3CH4  +  2O2  +  ION2,  pressurized  to  50-100  atm,  are  commoifly  used 
for  ram  accelerator  testing  at  ARL.  A  single  test  can  incorporate  several  different  gas  mixtures  in  different 
sections  of  the  accelerator  tube,  separated  by  diaphragms.  Fuel-rich  mixtures  are  used  with  fuel 
equivalence  ratios  as  high  as  3.  For  these  conditions,  a  good  understanding  of  CH4/O2  chemical  kinetics, 
especially  for  P  >  10  atm,  is  not  available  (Anderson  and  Kotlar  1991).  Accurate  numerical  simulation 
of  hydrocarbon-based  reacting  flow  is  very  demanding  in  terms  of  computational  resources  since  the 
number  of  intermediate  species  and  the  number  of  kinetic  steps  are  prohibitively  large  (i.e.,  ^15).  For 
high-pressure  systems,  the  computational  investment  may  not  be  justified  due  to  uncertainty  in  measured 
reaction  kinetics.  Thus,  global  reaction  mechanisms  (neglecting  intermediate  steps)  such  as 


CH4  +  2O2  CO  +  2H2  +  1.5O2 
CO  +  O.5O2  <=>  CO2 

2H2  +  O2  <=>  2H2O  (17) 

have  been  used.  The  reaction  rate  is  defined  using  the  Law  of  Mass  Action  and  an  Arrhenius  expression 
for  C,  the  specific  reaction  rate  constant. 


=  cnr.i  =  AT«exp 


a  b  c  d  e 
®CH4®02®C0®H2‘^C02  112 


^,0  ’ 


(18) 


7 


where  AT“  is  the  collision  frequency,  the  exponential  tenn  is  the  Boltzmann  factor,  is  the  activation 
energy,  and  Vj  represents  the  stoichiometric  coefficient  Exponents  a-f  are  chosen  such  that  the  reaction 
rates  match  results  from  flame  experiments  using  large  kinetic  mechanisms  (Westbrook  and  Diyer  1981) 
and  are  given  by  Nusca  (1993). 

An  Anhenius  reaction  rate  may  not  be  correct  for  the  high-pressure  flows,  and  much  of  the  necessaiy 
reaction  rate  data  for  hydrocarix)n  fuels  has  not  been  measured  above  10  atm  (Westbrook  and  Dryer  1981). 
In  addition,  the  three-step  kinetics  mechanism  described  previously  (Eq.  18)  is  based  on  the  low-pressure 
studies  of  Westbrook  (Westbrook  and  Dryer  1981)  where  the  exponents  are  determined  by  matching 
experimentally  measured  flame  speed.  For  the  ARL  full-bore  ram  accelerator,  where  mixture  ignition  is 
largely  determined  by  shock-induced  (i.e.,  reflected  shocks  from  the  accelerator  tube  wall)  and  viscous 
heating,  the  exponents  a-f  were  altered  to  the  degree  necessary  to  match  accelerator  tube  wall  pressures 
measured  at  ARL  (Nusca  1991,  1993,  1994;  Kraczynksi  and  Nusca  1992).  It  was  found  that  ram 
accelerator  performance  was  determined  to  a  laige  degree  by  the  reaction  induction  length  (i.e.,  ignition 
delay  time),  which  is  a  function  of  scale  and  mixture  (as  well  as  pressure).  For  the  preliminary  (TFD 
computations  described  in  the  present  report,  no  attempt  has  been  made  to  alter  die  kinetic  steps  or  the 
reaction  rate  expression  for  the  subbore  projectile  that  uses  ignititm  via  a  step  on  the  projectile. 
Experimental  confirmation  of  mixture  ignition  (i.e.,  location,  extent,  and  resulting  tube  wall  pressure 
history)  is  required  to  validate  the  kinetics  model  used  in  the  CFD  code.  However,  the  code  has  been 
extensively  validated  for  high-speed  compressible  gas  dynamics  (Palaniswamy,  Ota,  and  Chakravarthy 
1991). 

An  additional  issue  that  must  be  resolved  is  turbulence  modeling.  No  attempt  has  been  made  in  the 
presait  woik  (or  that  described  by  Nusca  [1994]  for  the  ARL  projectile)  to  test  the  dependence  of  results 
on  various  models  of  turbulence.  Nusca  (1993)  has  documented  flie  effects  of  viscosity  for  the  full-bore 
projectile  and  has  achieved  good  results  using  algebraic  turbulence  models.  No  attempt  has  been  made 
in  the  present  ^plications  to  alter  these  models. 

3.  RESULTS  AND  DISCUSSION 

Figure  2  shows  the  computational  mesh  for  ihe  projectile  with  a  step  located  at  the  forebody-midbody 
junction,  1  mm  in  height  The  forebody  is  conical,  85.1  mm  in  length  with  cone  half-angle  10°.  The 
midbody  is  cylindrical,  30  mm  in  length,  and  32  mm  in  diameter.  The  afterbody  is  conical,  98.6  nun  in 
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length,  and  with  a  half-an^e  of  3.5*^.  The  grid  was  generated  in  zones  so  that  sharp  body  junctions  are 
accurately  represented.  Approximately  11,000  grid  cells  cover  the  projectile  and  600  grid  cells  reside  in 
the  wake  (for  the  half-plane). 


A  gas  mixture  of  3CH4  +  2O2  +  ION2  was  used  for  the  CFD  simulations  described  in  this  section. 
The  mixture  is  characterized  by  a  frozen  speed  of  sound  a  =  361  m/s,  ratio  of  qrecific  heats  y  =  1.376, 
and  a  ChaiMnan-Jouget  detonation  velocity  V^j  =  1,450  m/s.  The  initial  pressure  and  temperature  of  the 
mixture  (before  projectile  injection)  was  50  atm  and  300  K,  respectively. 


Figure  3  shows  the  computed  frozen  flow  Mach  number  contours  in  grey-scale  (bright  white  represents 
zero  value,  i.e.,  the  projectile,  while  daik  black  represents  a  value  of  about  6)  for  the  iHojectile  geometry 
displayed  in  Figure  2  (step  height  of  1  mm)  and  a  fteestream  Mach  number  of  6  (V  =  2,175  m/s,  V/V^j  = 
1.5).  This  calculation  (as  are  all  calculations  described  in  this  report)  was  performed  a.s.suming  steady-state 
conditions.  Important  gas  dynamic  features  are  labeled.  A  bow  shock  ext^ids  from  the  projectile  nose, 
while  the  step  produces  a  normal  shock  and  a  small  stagnation  region.  The  merging  of  these  shocks 
produces  a  curved  "transmitted  shock"  (Figure  1)  that  reflects  from  the  accelerator  tube  wall,  downstream 
of  the  projectile.  Behind  the  transmitted  shock,  a  contact  surface  is  fonned  (not  iiiUy  visible  in  Mach 
number  contours)  that  separates  flows  of  different  density,  temperature  and  velocity  (i.e.,  flow  having 
passed  through  different  shock  waves).  The  recirculation  zone  behind  the  projectile  (wake)  is  clearly 
subsonic.  The  results  for  a  freestream  Mach  number  of  10  (V  =  3,610  m/s,  V/V^j  =  2.5),  displayed  in 
Figure  4,  show  a  downstream  bending  of  the  bow  and  transmitted  shocks  as  well  as  a  strengthening 
(i.e.,  producing  lower  values  of  the  downstream  Mach  number  flow,  greyer  color)  of  the  transmitted  shock 
and  step-generated  normal  shock. 


Continues  to  3  Projectile  Lengths  (J=1 95) 


Figure  2.  Computational  mesh.  Projectile  with  1-mm  sten. 
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Velocity  =  2175  m/s,  Mach  =  6.0,  V/Vcj  =  1.5 
Mixture:  3CH4  +  202  +  10N2,  a  =  361  m/s,  Gamma  =  1.376,  Vcj  =  1450  m/s 

Fill  Conditions:  P  =  50  atm,  T  =  300  K 
1  mm  Step,  3.5  deg  half-angle  afterbody 


Normal  Shock  Transmitted  Shock  Contact  Surface 


Tube  Wall  {D=120mm)  Reflected  Shock 


Figure  3.  Nonieacting  Mach  number  contours,  freestream  Mach  no.  6,  1-mm  step. 


Mach 


Velocity  =  3160  m/s,  Mach  =  10.0,  V/Vcj  =  2.5 
Mixture:  3CH4  +  202  +  10N2,  a  =  361  m/s,  Gamma  =  1.376,  Vcj  =  1450  m/s 

Fill  Conditions:  P  =  50  atm,  T  =  300  K 
1  mm  Step,  3.5  deg  half-angle  afterbody 
Normal  Shock  Transmitted  Shock  Contact  Surface 


Bow 

Shock 


9.658 

9.271 

8.885 

8.499 

8.112 

7.726 

7.340 

6.953 

6.567 

6.181 

5.795 

5.408 

5.022 

4.636 

4.249 

3.863 

3.477 

3.090 

2.704 

2.318 

1.932 

1.545 

1.159 

0.773 

0.386 


Figure  4.  Nonieacting  Mach  number  contours,  freestream  Mach  no.  10,  1-mm  step. 
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Figures  5,  6,  and  7  show  water  (H2O  is  a  major  product  of  the  combustion)  mass  firaction  contours 
in  grey-scale  (bright  white  representing  the  absence  of  water  and  dark  blade  representing  a  mass  fraction 
of  0.07  and  above).  These  calculations  are  an  extension  of  the  steady-state  results  drown  in  Figures  3 
and  4  with  the  flow  now  allowed  to  chemically  react  Figures  6  and  7  are  closer  views  of  the  step  region. 
Major  combustion  activity  is  predicted,  using  the  current  kinetics  mechanism  (Eq.  17),  in  the 
step-generated  stagnation  region  (behind  the  normal  shock)  and  continuing  around  the  step  and  then  over 
the  ix’ojectile  surface.  This  region  of  strong  combustion  and  combustion  products  effectively  resides 
between  the  omtact  surface  and  the  projectile  surface.  The  projectile  wake  flow  entrains  combustion 
products.  Note  in  Figure  6  that  the  production  of  HjO  is  reduced  somewhat  (i.e.,  lighter  grey  colors)  after 
expansion  around  the  step,  but  is  not  extinguished.  Figure  6  also  shows  a  close  view  of  the  same  step 
height  but  at  a  freestream  Mach  number  of  10.  The  H2O  mass  fraction  distribution  is  more  uniform  and 
extends  farther  into  the  flowfield  radially  away  from  the  projectile. 

The  bow  shock  wave  from  the  nose  does  not  initiate  dhemical  reaction.  The  detached  normal  shock 
wave  ahead  of  the  step  initiates  a  very  strong  detonation  wave  that  intersects  the  bow  shock.  At  the 
intersection  point,  a  transmitted  shock  is  directed  towards  the  tube  wall.  The  strength  of  this  transmitted 
shock  depends  on  the  flow  Mach  number  and  step  height.  The  transmitted  shock  weakens  as  its  distance 
from  the  projectile  increases,  imtil  reaching  the  tube  wall  and  reflecting.  At  the  intersection  point  between 
the  bow  shock  and  the  detonation  wave,  ahead  of  the  step,  a  contact-combustion  surface  is  also  generated. 
This  contact  surface  begins  at  the  intersection  point,  which  is  a  few  step  heights  above  ftie  surface,  and 
is  propagated  ftrom  there  t^atly  parallel  to  the  projectile  surface  towards  the  projectile  base  and  further 
downstream,  above  the  wake.  The  temperature  of  die  flow  with  combustion  products,  behind  the  stqi, 
is  about  5  to  8  times  the  freestream  temperature,  and  the  pressure  is  about  10  to  15  times  the  freestream 
pressure.  This  contact  surface  separates  the  hot  flow  of  the  combustion  products,  processed  by  the 
detonation  ahead  of  the  step,  and  the  outer  flow  that  passed  through  the  transmitted  shock  wave. 
Therefore,  the  contact  surface  acts  as  an  aerodynamic  cowling,  while  the  region  between  the  contact 
surface  and  the  projectile  body  surface  acts  in  a  manner  similar  to  a  ccxnbustion  chamber  in  a  conventional 
ramjet  engine.  However,  the  "combustion  chamber"  is  confined  by  the  aerodynamic  contact  surface 
without  requiring  the  stractural  cowling  of  a  conventional  engine.  The  hi^-ptessurc,  high-temperature 
combustion  products  in  this  aerodynamic  "combustion  chamber"  can  then  be  expanded  over  die  projectile 
afterbody  to  produce  thmst 
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H2O  Mass 

Velocity  =  2175  m/s,  Mach  =  6-0,  V/Vcj  =  1.5 
Mixture:  3CH4  +  202  +  10N2,  a  =  361  m/s,  Gamma  =  1.376,  Vcj=  1450  m/s 

Fill  Conditions:  P  =  50  atm,  T  =  300  K 
1  mm  Step,  3.5  deg  half-angle  afterbody 


Tube  Wall  (D=120mm) 


m 


Fraction 

0.070 
0.067 
0.064 
0.061 
0.058 
0.055 
0.052 
0.050 
0.047 
0.044 
0.041 
0.038 
0.035 
0.032 
0.029 
0.026 
0.023 
0.020 
0.018 
0.015 
0.012 
0.009 
0.006 
0.003 
0.000 


Figure  5.  Reacting  H2O  contours,  freestream  Mach  no.  6,  1-mm  step. 


Figure  6.  Reacting  H2O  contours,  freestream  Mach  nos.  6  and  10,  1-mm  step. 


12 


Reacting  Flow  Field 

Velocity  s:  2175  mis,  Mach  s  6.0,  V/V<^  =  1.5 


HjO  Mass  Fraction 

““  0.070 


Reacting  Flow  Field 

Veloe% = 2175  ntfs,  Mach = 6.0,  V/Vq= 


H3O  Mass  Fraction 

0.^70 


Figure  7.  Reacting  HoO  contours,  frecstream  Mach  no.  6,  0.5-  and  1.5-mm  steps. 


Figure  7  shows  water  mass  fraction  contours  in  the  vicinity  of  the  0.5-  and  1.5-mm  steps  for  a 
fieestieam  Mach  number  of  6.  The  H2O  mass  fiaction  distributions  are  similar  to  those  for  the  1-mm  step 
height,  but  with  the  combustion  region  remaining  close  to  the  projectile  body  for  the  0.5-mm  step  and 
extending  farther  into  the  flowfield  for  the  1.5-mm  step.  The  larger  step  causes  a  more  extensive 
stagnation  region  in  front  of  the  step  (behind  the  normal  shock). 


As  the  step  hei^t  is  increased,  the  conditions  and  the  position  of  the  contact  surface  are  effected, 
resulting  in  increased  mass  of  combustion  products  at  higher  pressures  and  temperatures,  but  with  an 
increase  in  drag  due  the  increased  step  height.  The  observation  that  combustion  can  be  initiated  and 
sustained  for  the  0.5-mm  step  is  in  general  agreement  with  the  approximate  analysis  of  Rom  and  Avital 
(1992)  where  the  minimum  step  height  for  establishing  a  combustion-detonation  front  ahead  of  a  step  for 
stoichiometric  CH4/O2  mixture  is  evaluated.  This  step  height  is  about  0.7  mm  for  Mach  no.  6  and 
precombustion  pressure  of  1  atmosphere. 

The  computed  pressure  distributions  on  the  projectile  surface  were  integrated  to  form  a  thrust 
coefficient,  F/PA,  where  P  is  the  initial  mixture  pressure  (50  atm)  and  A  is  the  cross-sectional  area  of  the 
accelerator  tube  (diameter  of  120  mm).  For  a  freestream  Mach  number  of  6  and  the  1-mm  step,  the 
nonreacting  flow  yielded  F/PA  =  -0.4  (i.e.,  drag),  while  the  reacting  flow  yielded  F/PA  =  3.1.  Projectiles 
with  step  heights  of  0.5  mm  and  1.5  mm,  and  nonreacting  flow  yielded  F/PA  =  -0.3  and  -0.9, 
respectively.  With  chemical  reactions,  these  cases  yielded  F/PA  =  3.4  and  3.06,  respectively.  The  0.5-mm 
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step  yields  a  smaller  drag  that  must  be  canceled  by  the  thrust  produced  by  combustion  and  thus  yields  a 
larger  thmst  coefficient  The  1.5-mm  step  yields  a  larger  drag  that  must  be  canceled  by  the  thrust 
produced  by  combustion.  Nevertheless,  the  more  extensive  combustion  for  this  step  height  yields  a 
positive  dmist 

The  author  feels  confident  about  the  fidelity  of  the  CFD  code  with  regards  to  gas  dynamic  ph^omena 
displayed  in  this  section.  However,  the  accuracy  of  the  assumed  chemical  kinetics  mechanism  can  be 
argued,  especially  at  high-pressure  and  fuel-rich  conditions.  Further  research  and  validation  is  warranted. 
The  presented  results  should  therefore  be  considered  preliminary. 

4.  CONCLUSIONS 

The  present  numerical  simulation  of  the  effects  of  a  forward-facing  step  positioned  on  the  body  of  a 
hypervelocity  projectile,  clearly  indicates  that  the  step  wiU  produce  chemical  reactions  of  the  hydrocarbon 
atmosphere.  The  present  calculations  also  show  that  these  reactions  do  not  extinguish  after  flow  expansion 
around  the  step,  ate  sustained  over  the  afterbody  and  the  wake  of  the  projectile,  and  can  be  used  to 
produce  thrust.  It  is  ^own  that  the  interaction  between  the  bow  shock  from  the  ccmical  fotebody  and  the 
detached  shock-detonation  wave  from  die  forward-facing  step,  which  is  positioned  on  the  shoulder  of  the 
projectile,  results  in  confinement  of  a  combustion  layer  about  the  projectile  afterbody.  The  high-pressure 
and  temperature  products  of  the  combustion,  which  are  confined  to  this  layer,  are  then  eiqianded  into  the 
rear  and  the  base  flow  of  the  projectile  and  can  be  used  to  generate  thrust.  However,  since  die  projectile 
sh^  and  combustion  scales  are  different  in  the  present  external  propulsion  case  from  those  in  the  Ram 
Accelerator  case,  these  results  should  be  validated  by  an  experimental  firing  of  the  projectile  geometry 
investigated. 
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LIST  OF  SYMBOLS 


A  cross-sectional  area  of  launch  tube 
b  covolume 

Cp  specific  heat  capacity,  constant  p 

Cy  specific  heat  capacity,  constant  volume 

C  specific  reaction  rate  constant 

D  mass  diffusion  coefficient 

e  specific  total  internal  aiergy 

activation  energy 
F,G  flux  vectors  (Eq.  1) 

g  Gibbs  energy  per  mole 

h  molar  specific  enthalpy 

i  i-th  species 

L  total  body  length 

m  mixture  quantity 

M  molecular  weight 

n  number  of  moles 

N  total  number  of  species 

p  static  pressure 

R  specific  gas  constant,  (y-l)Cp/V 

universal  gas  constant,  RMjj, 

Re  Reynolds  number,  pUL^i 

Sc  Schmidt  number,  Pn/pD 

t  time 

T  static  temperature 

u  axial  velocity 

V  radial  velocity 

U  magnitude  of  the  local  velocity  vector 

V  specific  volume  (1/p) 

W  dependent  variable  vector  (Eq.  1) 

x,y  cartesian  coordinates 

X  species  mole  fraction 
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a  0  for  2D,  1  for  axisymmetric 

Y  ratio  of  specific  heats, 

AHf  enthalpy  of  formation 

K  heat  transfer  coefficient 

p  molecular  viscosity 

V  stoichiometric  coefficient 

p  density 

o  species  mass  fiaction 

T  shear  stress  tensor 

<D  chemical  production  term  (Eq.  18) 

source  term  vector  (Eq.  1) 
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